A Kaluza-Klein Subttractor 

Sanjib JANAfjand Chethan KRISHNANf] 



a Center for High Energy Physics 
Indian Institute of Science, Bangalore, India 



Abstract 



We generalize the results of arXiv:1212.1875 and arXiv:1212.6919 on attraction basins and 



their boundaries to the case of a specific class of rotating black holes, namely the ergo- 
free branch of extremal black holes in Kaluza-Klein theory. We find that exact solutions 
that span the attraction basin can be found even in the rotating case by appealing to certain 
symmetries of the equations of motion. They are characterized by two asymptotic parameters 
that generalize those of the non-rotating case, and the boundaries of the basin are spinning 
versions of the (generalized) subttractor geometry. We also give examples to illustrate that 
the shape of the attraction basin can drastically change depending on the theory. 
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1 Introduction 



Attractor mechanism pQ is the statement that the moduli scalars in an extremal black 
hole geometry approach fixed values at the horizon even though their radial profiles can be 
different. In some recent work [21 E] we demonstrated that the boundary of the basin of 
attraction of a class of static attractor black holes in dilatonic supergravity is the so-called 
(generalized) subtracted geometry [H |5j [6] in the extremal limit. We called these solutions 
subttractoraS 

In this paper, we wish to generalize this statement to a class of examples where the 
extremal black hole is rotating, and determine the boundary of the attraction basin in this 
more general CclSG clS well. This is non-trivial because adding rotation breaks the spherical 
symmetry of the problem down to an axial symmetry. 

The theory we will work with is the simplest Kaluza-Klein theory, namely Einstein- 
Maxwell-dilaton theory obtained by reduction of the 5D Einstein action [TU]. In this paper 
we will study (rotating) extremal dyonic black holes of this theory with electric and magnetic 
charges. In section 2, we will show how in the static case, the equations of motion of this 
theory can be fully integrated and the extremal solutions are characterized by two parameters 
(integration constants) which we call d\ and d 2 . The basin of attraction of asymptotically 
flat dyons is described by a certain domain in the di,d 2 space (which we describe). This 
domain of d\,d 2 that we find here is quite different from the attraction basin found in the 
theories in [21 [3]. To emphasize that the domains of attraction and the structure of the 
solutions can vary drastically, in an appendix)^] we also discuss this domain for yet another 
fully solvable system (case III in [TTj ) and describe its basin structure. 

In section 3, we will consider the same Kaluza-Klein theory, but now we will look at what 
we are really after: rotating (extremal) solutions. Once we turn on rotation, the system is 
substantially more complicated. The black hole solutions of this theory were written down in 
[T2| [131 E]. There are two kinds of extremal limits that these black holes have, the so called 
ergo-branch and the ergo-free branch [15J. The branch that allows a natural generalization 
from the static case is the ergo-free branch, and that is the one we will study. As it turns 
out, the black hole solutions of [121 [131 EH] are not general enough to discuss the attraction 
basin that generalizes the static basin of section 2. Fortunately, despite the complications of 
the rotating case, using some symmetries of the equations of motion we are able to generate 
precisely such solutions. These new solutions are again characterized by two parameters 

1 Subtracted geometries initially arose in the context of the Kerr/CFT correspondence [7] and as a geo- 
metrical way to manifest the hidden conformal symmetry of black holes [5] . Subtracted geometry has been 
constructed and further discussed in [JJ [SJ [HI [21 131 M ■ 

2 The theory considered in this appendix is a dilatonic toy model, independent of the KK theory we discuss 
in the rest of the paper. 
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which we call di, d2 in analogy with the static case. 

We conclude in section 4 by discussing the basin of attraction of these new solutions 
and presenting the full attractor flow structure for various values of 9, the polar angle in 
the geometry. The boundaries of attraction are characteried again by certain limits of the 
parameters d\ and di. The discussion of these boundaries is somewhat distracting, so we 
have relegated it to an appendix. In these limits, the solutions degenerate and turn into a 
generalization of the (rotating) subttractor geometry. The warp factor of the geometry goes 
as ~ r 3 for these solutions, which should be contrasted to ~ r 4 for flat space and ~ r in the 
original subtracted geometry of jU [6] . 

Note added: After this work was substantially completed, two papers appeared on 
the arXiv which deal with related ideas. [16] constructed^] general black hole flow solutions 
and [IB] deals with extremal black holes in dilatonic supergravity. These works rely on the 
Harrison transfromation machiner j^j to generate solutions. These solutions should enable an 
exhaustive study of the attraction basin structure of general (extremal) black holes [21] in 
these theories, going far beyond our approach here, which is based on symmetry arguments 
and other parlor tricks. 

2 The Static Attractor: Kaluza-Klein Dyon 

The reduction of 5D Einstein action on a circle gives [10] 

S = J d 4 xy^(R-2(d(j>) 2 -e^F^F^. (2.1) 

In this section we will look for general static extremal solutions of this theory. The static 
"attractor ansatz" that we will use is 

ds 2 = -airfdt 2 + ^— + 6(r) W, (2.2) 
a{r) 2 

e -2 v / 30 

F = Qi sin 6d6 A d<f) -\ — — Q 2 dt A dr, <p = <p(r). (2.3) 



The resulting equations of motion are 



(a 2 b 2 )" -2 = (2.4) 
J 



'"" +0' 2 = O (2.5) 



,.2,,/y ^Kff(0) 



(aW/ - = (2.6) 



3 Generalizing [17]. 

4 Scc [19, 20, 9J for previous discussions on Harrison transformations in the subtracted geometry context. 
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and the energy constraint from Einstein equations is 




a 



2 



+ 



VM) 
b 2 



- a 2 b 2 4>' 2 -1=0. 



(2.7) 



The effective potential V e g is 



Kff(0) = e 2 ^Q? + e 



2 v / 3</>q2 



(2.8) 



Here, we obtained the equations of motion and the effective potential by using an electric 
plus magnetic (ie., dyonic) ansatz in a theory with one gauge field. It is straightforward to 
see that one can also get these exact same equations, by (say) starting with a theory where 
instead we have two gauge fields with suitable couplings to one scalar, with only magnetic 
charges turned on. 

These equations can in fact be solved exactly [TTI |2"2"| |2"3"] . This has been done in the 
appendix of pj] after imposing the condition that the system has a double-zero horizon (ie., 
it is extremal). But they make the further assumption that the asymptotic value of a(r) is 
unity, which is how asymptotic flatness is imposed. We find it useful for our purposes to 
leave this freedom arbitrary. This has a few virtues. One is that this essentially gives us the 
most general solution of the system which is extremal: in other words all the near-horizon 
data has been fixed, but none of the boundary data is. The second virtue is that in many 
cases, we need this degree of freedom to gain a full understanding of the attraction basin. 
Finally, it is evident from the equations of motion above that there is a rescaling symmetry 
to the system under 



In effect, fixing the value of a(r) at r — > oo fixes this scale. Not fixing this scaling is 
tantamount to leaving it as an integration constant in the solutions. This is what we will 
do, because a generalization of this scaling symmetry exists also in the rotating case as well 
and we will make use of it in generating our new solutions there. This is crucial because 
unlike in the static case, solving the system frontally is an industry in itself in the rotating 
case: the equations are complicated partical differential equations and the system has less 
symmetry. 

It is straightforward to construct the general solution by adapting the results in Appendix 
B.3.2 of [TT]. The result can be expressed as 



r — > r/fj,, a — > a/ ' \x, 



(2.9) 



i^ = 9l l 2(d 2 -d 1 )r 2 -2d 2 r + l 
~~ Qi [ 2d 1 r 2 - 2d 2 r + 1 



(2.10) 




2QiQ 2 v/(2(d| - d\)r 2 - 2d 2 r + l)(2dir 2 - 2d 2 r + 1) ' 
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Figure 1: The attraction basin of the Kaluza-Klein black hole. The basin structure in terms 
of d\ and d 2 is identical, whether the black hole is spinnning or not. The hatched region is 
bounded by the curves d\ = di and di = 0. 

Here d 2 and d\ are the integration constants, and the asymptotic values of a(r) and </>(r) 
(lets call them do and (p^) can be immediately read off from the above expressions in terms 
of them. The results of [TT] were obtained for the special case, oq = 1. 

Now, the attraction basin structure can be easily studied using the approach spelled out 
in [3]. The idea is that we want the solutions to be regular in the region < r < oo and 
this will happen only if the quadratic polynomial expressions that show up in the solution 
above have no zeros in that range. The zeros (in r) of the polynomials can be determined in 
terms of di and d 2 and this imposes inequalities on the ranges of values that d 2 and d\ can 
take. We present the analyses of these inequalities in an appendix, the final result is that 
the attraction basin is the acute region 

di > 0, d 2 < 0, and d\ > d x , (2-12) 

in the d\ — d 2 plane that is bounded by the curves 

d l = 0, di = dj, with d 2 < 0. (2.13) 

This is the hatched region in figure 1. Note that the attraction basin is no longer a quadrant 
as it was in the cases considered in [3 J . To emphasize that the curves defining the attraction 
boundaries can take other forms, we consider another exactly solvable dilatonic gravity 
theory in an appendix and describe its attraction basin in the static case. 

If we slice through the attraction basinQ along a d 2 = constant line, we will see the various 

5 We will use the phrase "attraction basin" to refer to the hatched region in the d\ — di parameter space. 
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asymptotically flat hairy solutions that arise as we go from one boundary of the attraction 
basin to the other. The resulting plot is identical to the plots with /3 = 0or# = 7r/2in figure 
2, so we will not repeat them here. If we perturb a black hole with parameters d\ = d\ and 
d 2 = d 2 near the horizon by spherically symmetric perturbations and evolve it radially, we 
will find another solution in the same family with a different value of d\ and d 2 . In fact this 
new solution and old solution will lie precisely on the line 



d 2 = negative constant 



4 



(2.14) 



in the attraction basin. The perturbation theory is straightforward (even if tedious) to do by 
a simple adaptation of the results in Appendix D of [3] (see also very similar computations 
in [2U |25l [26] for some more detail). We will not emphasize it here, because we are merely 
interested in some illustrative slicing of the attraction basing The only salient point is that 
because of the specific form of the couplings in our system, the leading non-zero term in 
in the scalar perturbation appears not at 0(r), but 0(r 2 ), see eg. [2j. It turns out that 
this perturbation (lets call it 02 ) is the only degree of freedom for radial fluctuations, and it 
therefore parametrizes a line along the basin: 



di = d\ 



<P2 



do = <£. 



(Trivially) eliminating 02 from system is how we obtained (2.14). 

At the boundaries of attraction, the solutions degenerate. At the d\ 
becomes 



(2.15) 



broundary it 



02 

Qi 



2d\r 2 



2d 2 r + 1 



a 2 (r) 



2d 2 r 



3/2 



2QiQ 2 ^{2d\r 2 - 2d 2 r + 1)(1 - 2d 2 r) 



b(r) = r/a(r) 



(2.16) 
(2.17) 



and at the d\ 



d 1 end, the solution is 



Qi 



1 - 2d 2 r 



2d\r 2 



2d 2 r + 1 



3/2 



a 2 (r) 



1 



2QxQ 2 - 2d 2 r)(2d 2 2 r 2 - 2d 2 r + 1) 



b(r) = r/a(r) 



(2.18) 
(2.19) 



The warp factor A defined by v^A = r 2 /a 2 determines the asymptotic behavior of the 



geometry, and it goes as 



in both cases. So the solutions in this limit stop being 



asymptotically flat, and are examples of subtracted geometries [H El [13, [2j [3] . Note that the 
two solutions above are essentially identical except for a sign flip (and shift) in the scalar. 



6 The slicing emerging from the perturbation theory is perhaps the most canonical slicing of the attractor 
basin, however. 
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3 Adding Rotation: the Ergo-Free Branch 



Adding rotation makes the black hole significantly more complicated. Rotating Kaluza- 
Klein black holes were first constructed in [121 E31 EI]- In the non-extremal case, on top of 
a mass parameter and the charges Q\ and Q2, now we have a rotation parameter a& in 
the metric \IE\- Since our interest is in the attractor mechanism, we will only be interested 
in extremal solutions and will not write down the general non-extremal metric. 

Starting with the non-extremal case, there are two kinds of extremal limits one can 
take. One involves setting = Mk- This is called the ergo branch because the solution 
has an ergosphere. The other extremal limit is taken by sending ajt, — y while holding 
dk/Mk = (3 fixed. This limit is called the ergo-free branch and this is what we will be 
interested irj^J Regularity of the solution dictates that \/3\ < 1 [12] • Note that the sign of /3 
merely captures the direction of spin. 

Our starting point will be the hairless extremal solutions discussed in [121 EH E2] . We 
will write the metric using the electric and magnetic charges (Q, P) in the normalization of 
[12] . In the ergo-free extremal limit, and Mk do not show up because they have gone to 
zero, but their ratio (3 doees. The metric takes the form: 



ds 2 - ~^=(dt - ujd^f + ^P^r 2 + ^/fj 2 dn 2 (3.1) 
v/1/2 A 



with 



2 P*/*y/PW + QW P^Q^ 

fi=r + — p= r + (l-0cos0), (3.2) 



2 Q2/3 v /p 2 /3 + g2/3^ p2/3g4/3 ( 



/ 2 = r 2 + ^ v - ^ r+ ~ (l + ffcosfl), (3.3) 

'7T S7T 



2 /3sin 2 g 

X = r , w = . 3.4 

87r r 

The lack of spherical symmetry entrers the metric through the non-trivial angular depen- 
dence on the time fibration. 

In the static case we could use an ansatz for the field strength which automatically solved 
the gauge field equations without using the details of the specific solution (cf. the attractor 
ansatz from the previous section). But now, we have to work with gauge fields directly and 



7 In the non-rotating case, the extremal limit in these coordinates corresponds to sending Mk — y 0. So 
the ergo- free branch is the natural generalization of the static extremal limit. Understanding the attraction 
basins of all possible extremal solutions is likely to require full control on all the hairy solutions. This should 
be possible to do in light of the recent results in [TBI UH] , but we will not attempt it here. 
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the non-zero components take the following form: 

Q 



4v^ h 



p2/3^p2/3 + g 2 /3 

r+ ± - — l + /3cosi 



P 



Together with the scalar 



cos 9 + 



P Sin 2 9 f Q 2/3 v /p2/3 + g2/3^ ^ p2/3g4/3 ^ 



2/ 2 V 



47T 7 



exp 



(2^30) 



(3.5) 



(3.6) 



(3.7) 



this completes the description of the black hole. 

Note that the solution above has no free parameters other than the charges of the hole 
(counting also (5). What we would like is a generalization of this solution which has two 
free parameters (integration constants) generalizing what we saw in the static case, so that 
we can conveniently describe the attraction basin. One of these parameters corresponds to 
the freedom to scale the asymptotics of the metric and the other one is the freedom to shift 
the asymptotic value of the scalar. The latter is simple to implement [TS]. We see from the 



action (2.1) that the theory has a symmetry under 

+ A 



(3.8) 



This transformation however changes the charges of the black hole, but we can rescale the 
charge parameters appearing in the solution while doing the above shift /scaling to obtain a 
one-parameter family of solutions with identical electric magnetic charges |15j : 



+ A, 



-V3A 



(3.9) 



This is a (one parameter) solution generating transformation that we can use. 

But where does the second parameter come from? The hint here is that the static ansatz 
had a scaling symmetry in a which could be removed by rescaling r. We will look for a 
generalization of this to the rotating case here. But since we do not have an ansatz when 
rotation is present^ we cannot immediately search for this symmetry in the equations of 



motion. But an inspection of the metric (3.1) and the gauge field 



A = A t dt + A v dip 



(3.10) 



8 Naively adding arbitrary 9 dependence to the ansatz functions results in PDEs. Moreover, choosing an 
ansatz for the gauge field will be non-trivial. However, even though we have not explored this in detail, it 
seems possible that by choosing the 9 dependence of the ansatz carefully and using inspired analogies from 
the exact solution EH [33] , one might get a tractable set of EoMs. The solution we find here is enough 
for our purposes since it is applicable in the extremal limit, but the static "attractor" ansatz is applicable 
even away from extremality. 
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reveals that they (and trivially, the scalar) are invariant under 

r ->■ r/fi, t^tfi, X ->■ X//J 2 , u ->■ cu /x, A t ->■ (3.11) 

This means that if we have a solution, replacing r by r//i will result in a new solution if we 
compensate for it by an opposite scaling of u,X and A t . 

Doing the above two transformations results in the new metric functions: 



h = r 2 /f + 
h = r 2 /^ 2 + 



P^P^ + QT 


r/fi + 


p 4/3 n 2/3 




8tt 




r/pL + 


p 2/3 n 4/3 

-*o V o 




8tt 


X = r 2 , w = 


PoQo 
8tt 


/3 sin 2 6 
r 



(1 -/9cos0), 
(1 + /3 cos 9), 



(3.12) 



(3.13) 



(3.14) 



The gauge field components (in terms of the new metric functions fi and f 2 ) take the form: 



A t = -/I 



4v^ / 2 



r//i + 



o + Qo 



2/3 



4^ 



1 + /3cos( 



A, 



P e 



-V3\ 



4A 



cos + — — — I - j= r/n + 



2v^F 



Aix 



(3.15) 



(3.16) 



Finally, the scalar (again in terms of the new f\ and f 2 ) is: 

exp( 2^30) = exp(2^3A) 



Ml 
h 



(3.17) 



It is important to note that in the above expressions, P$ and Qo are the scaled charges: 

Qo = e-^ x Q, P = e^ x P. (3.18) 

These are the two parameter solutions we seek, and they are pretty ugly and complicated in 
terms of (3, P, Q, \i and A. 

But since we have two free parameters at this stage, we can compare their static limit 
with the static solutions of section 2 and write them in a cleaner notation. It turns out that 
this notation also lends itself to a better description of the attraction basin. Comparing the 
expressions after setting (3 = in the formulas above, results in the relations 



A 



1 



2v/3 



log 



Q 2 (d\-d 1 \^ 



Qi v di 



n l/4 



di(d\ - di) 



(3.19) 



P = -VttQi, Q = -4v^Q 2 . 



(3.20) 
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So apart from a minor difference in the normalization, the charges P, Q of the black hole 
solutions of [121 [15] are the same as in our notation, Qi, Q 2 - In terms of the (Qi, Q 2 , di, d 2 ) 
notation, then the full rotating solutions take the much more tractable form 



X(r) = r , u> 



2 2 Qi Q2 o . 2 /j 

■ — psm 0, 

r 



/1 = 2 Qi Q 2 



/ 2 = 2 Qx Q s 



G?o — d\ 



d\ 



{2d 1 r 2 -2d 2 r + l-f3 cos9), 



d x 



di — d\ 



{2[d\ - rfi)r 2 - 2d 2 r + 1 + cos 0) 



with the gauge fields 



.4, 



Qi / <h 

2f 2 V(dl-d 1 



2d x r - d 2 (l + /3cos( 



2QxQ 2 / c?i . 2 

cos# — w ~2 — — (a 2 r — ljpsm v 

h V "2 ~ "1 



and the scalar 



exp(2v / 30) = ^(^^j 



(3.21) 
(3.22) 

(3.23) 



(3.24) 
(3.25) 



(3.26) 



We have checked explicitly that these comprise a solution of the fulr] set of equations of 



motion arising from the action (2.1). This is the form of the solution that we will use to 



investigae the attraction basin in the next section. 



4 The Basin 

The solutions presented in the previosu section are regular only for specific ranges of the 
parameters d\ and d 2 . These ranges can be detremined by systematic application of the 
requirement of regularity, and is done in an appendix. The end result that one finds is what 
we will call the basin of attraction. It turns out that the basin for the rotating black hole 



for any value of 9 is bounded by the same curves as in the static case (2.12 2.13). We plot 
the solutions for a few values of /3 and 9 in figure 2. As expected, the plots with 9 = n/2 or 
(3 = give identical curves - these correspond to the attraction basin of the static attractor. 

It is worthwhile making a couple of comments about these solutions. The attraction 
basin structure is determined by the regularity of the functions fx and f 2 (see Appendix A). 



3 As opposed to a solution of the attractor ansatz of section 2, which is only valid in the static limit. 



9 



Figure 2: Curves in each row is a fixed /3(= 0, 0.5, 1, left-right), and each column is a fixed 
9(= 0, 7r/4, 7r/2, top-down). The last non-divergent line is the attraction boundary. 



It is clear from the form of these functions that when 



(3 cos 9 — > —(3 cos 9, 



this is tantamount to the switch 



d\ — > d\ — di, d 2 — > d 2 . 



(4.1) 



(4.2) 



But under such a replacement, the attraction basin (2.12 2.13) retains its structure. This 



means that we only have to determine the attraction basin for positive values of f3 cos 9. 
At the boundaries of attraction: 



hh = 4QlQi(-2d 2 r + 1 -0cos9)(24 r 2 - 2d 2 r + 1 + (3 cos 9), 
Q 2 (1d\ r 2 -2d 2 r + l + /3cos^ 



exp 



(2\/3» 



Qi\ -2d 2 r + 1 - (3 cos 9 
when d\ — 0, and 

hh = ±Q\Q\{2d\ r 2 -2d 2 r + l-f3 cos 9)(-2d 2 r + 1 + (3 cos 9), 



exp 



(2^3>) 



Q 2 / -2c? 2 r + 1 + g cos 9 
~Q[\2d 2 2 r 2 -2d 2 r + 1 - (3cos9 
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(4.3) 
(4.4) 

(4.5) 
(4.6) 



when d\ = d\. From these it follows that the metric and the scalar both have well-defined 
limits at the boundaries. Asymptotically, the warp factor goes ~ r 3 as in the non-rotating 
case, so this is a rotating version of the (generalized) subttractor geometry found in [21 E]- It 
is worth noting here that the solutions here again exhibit a symmetry under the flipping of 
the two boundaries, essentially upro a flip in the sign of the scalar and (3 cos 9 —> — (3 cos 9. 
This is related to the discussion in the last paragraph. 

We haven't explicitly written the gauge fields, but they are easily computed. On the 
d\ = d\ boundary the A t diverges and this might cause some consternation. But this 
divergence is actually a gauge artifact: the malady arises from a divergent constant and one 
can explicitly check that the field strengths remain finite at both boundaries. 



Acknowledgments 

We thank Avik Chakraborty and Amitabh Virmani for discussions, and Sudhir Vempati 
for teaching us a simple trick that made it possible for Mathematica to do its job without 
choking. 



A The Boundary of the Basin 

Determining the boundary of the basin is in principle straightforward: one just has to 
demand the functions that define the fields should be regular and real for < r < oo. 
But the problem splits up into many cases, so it is worthwhile being systematic. We will 
determine the basin here for the case of the rotating solution we constructed in section 3, 
we can reproduce the results we stated in the static case of section 2 by demanding that (3, 
the rotation parameter, is zero. 

The first thing to note is that regularity is determined by the functions fx and f 2 defined 
in section 3. Regularity is equivalent to the statement that these functions should only have 
zeros away from the positive real axis. To avoid clutter, we will temperorily introduce the 
functions: 

f x = ai r 2 + h r + 1, (A.l) 
f 2 = a 2 r 2 + b 2 r + 1, (A.2) 

and phrase our discussions in terms of them. At the end of the discussion 10 we will relate 

10 We will occasionally suppress the index i in what follows, when we are making general statements valid 
for either i. 
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di=i,2 and 61=1,2 to d\ and d 2 . These relations can be taken to be 

ax = — , 61 = , (A.3) 

1 — a 1 — a 

a 2 = -^, = -—?-. A.4 
1 — a 1 — a 

where we have intrdoduced 

a = /3cos# (A.5) 

for convenience. Since /3 = ax/Mk < 1 and |cos#| < 1 we have the condition that —1 < 
a < 1. 

Before we proceed, we make one comment: the end result for the attraction basin is simple 



n 



and elegant, but embarassingly, we have not been able to derive it using a simple argument 
It seems evident that our result should follow from high-school notions about polynomials, 
roots, positivity, etc. but in the following, we have found it quicker to brute-force our way 
to the answer case-by-case. 

First we notice that both ai and a 2 have to have the same sign. If not, for some large 
enough value of r, both /1//2 an d f\f 2 will become negative and run into trouble with the 
square roots in the solution. 

Now, if a < 0, the roots of / cannot be complex: because D = b 2 — 4a will be forced to 
be non-negative. On the other hand, when roots are real, we want them both to be negative 
so that there are no zeros on the positive real axis. For a < 0, this means that both 

- b - Vb 2 - 4a > 0, and - b + Vb 2 - 4a > 0. (A.6) 

For negative a, it is clear that these equations cannot both be be true at the same time. 
This means that we must have a positive. 

Now we discuss the case a^ > in detail. There are two ways in which the function / 
can be regular on the entire positive real axis: it can either have compex (and mutually 
conjugate) roots, or both its roots are on the negative real line. Complex conjugate roots 
happen when D = b 2 — 4a < 0. On the other hand, if the equations have real roots (ie., 
D > 0), the bigger one will be 

1 



'-b + y/Di). (A.7) 



11 Strictly speaking, our case-by-case argument gets complicated only when there is rotation. When there 
is no rotation, many of the separate cases collapse onto each other. Interestingly, the final result for the 



basin in terms of the — d\ parameters is the same for both static and rotating black holes ( 2. 12||2. 13 1. In 
particular it does not depend on the polar angle 9 or the rotation parameter /3. 
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We want this root to be < 0. This forces b > 0. 

So in effect, what we have to do in order to determine the domain of attraction (the 
domain where the solutions are regular) is to consider the various possible combinations 
arising from f x and In principle there are four such cases (a« > has to be also imposed 
for all of these cases as we already dicussed): 

• b\ - 4ai < with b\ - 4a 2 < 0, 

• b\ - 4ai < with b\ - 4a 2 > 0, b 2 > 0, 

• b\ - 4ai > 0, bi > with b\ - 4a 2 < 0, 

• b\ - 4ai > 0, bi > with 6^ - 4a 2 > 0, b 2 > 0. 

We will call these the four branches of the basin. We will take the intersection of the regions 
in each bullet, and then take the union of the four regions to get the full attraction basin. 



In terms of (A. 3 A. 4), these four branches take the form 



Branch 1: 



^(2a-l) + 2di(l-a)<0, d\ - d x > 0, 
dl - 2di(l + a) < 0, d x > 



(A.8) 
(A.9) 



Branch 2: 



d 2 2 {2a - 1) + 2d 1 (l - a) < 0, d\ - d x > 0, 
da - 2di(l + a) > 0, d 2 < 0, d x > 0, 



(A.10) 
(A.11) 



Branch 3: 



d 2 2 (2a - 1) + 2di(l - a) > 0, d 2 < 0, d\- d x > 0, 
- 2di(l + a) < 0, di > 0, 



(A.12) 
(A.13) 



Branch 4: 



d^(2a - 1) + 2di(l - a) > 0, d 2 < 0, d\- d x > 0, 
d^ - 2di(l + a) > 0, d 2 < 0, di > 0, 



(A.14) 
(A.15) 
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Figure 3: The contribution to the basin from branch 4, for case 3. The curves that are 
plotted are the boundaries of the inequalities in the branch, and their intersection is the 
blackened region. More concretely, the branch 4 inequalities in terms of 7 = —a take the 
form g?2 < 2d\ yj^^ j , d\ > 2d%(l — 7), < 0, d\ > di, d\ > 0. We have made the plot above 
for the specific value 7 = 0.4. The ordering of the curves is fixed by the fact that < 7 < |. 



At this stage we have to distinguish four cases again, because the qualitative features of the 
various branches depend on the ranges of a\ These are 

Case 1 : < a < -, Case 2 : i < a < 1, (A. 16) 

Case 3 : - - < a < 0, Case 4 : - 1 < a < --. (A. 17) 

For each of these cases, one can investigate the four branches of the solution. The rest of the 
computation is merely an exercise in persistence, so we will merely present one illustrative 
subcase as an example and then quote the final result. The specific case we present is Branch 
4 of < 7 = —a < |. The contribution to the basin from that branch, arising from the 
intersection of the various inequalities, is presented in figure 3. Note that in the figure, the 
ordering of the curves depends crucially on the value of 7 and this is the reason why we had 
to work with the four cases separately for the ranges of a. 

Once the dust settles, in each four of the ranges of a, despite the vastly different inter- 
mediate steps, the union of all the four branches results in the same attraction basin: 

d 2 2 -d 1 > 0, d 2 < 0, di > 0. (A.18) 

In particular, this result is a- independent for all values of \a\ < 1. The boundaries of 
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attraction are obtained when one replaces the first two inequalities with equalities, and 
retains the last one as an inequality. 

In principle, we need not have computed the basin separetely for negative a, once we 
compute the basin for the range < a < 1. This is because of the symmetry presented in 
(4.1 4.2). However, the intermediate results in the computation for positive and negative a 
are (at least superficially) quite different, so this works as a check of our results. 



B One More Static Example 

In this section, we will explore a case where the equations of motion are of the form 



(2.4 2.7) but with the effective potential given by 

Kfr = e 40 Ql + e" 60 Q\. (B.l) 

Together with the resutls in |3j and the main body of this paper, this shows that the attrac- 
tion basin can change completely depending on the system we consider. 

The above specific form of the effective potential is another example where the resulting 
EOMs are integrable. This effective potential arises, for example, in an Einstein-Maxwell- 
dilaton system with one scalar and two vector fields (with magnetic charges Q\ and Q 2 
turned on for the gauge fields) , with a gauge coupling matrix (see [ITJ [2j [3] ) given by 



fM =[ n _™ • (B.2) 




By relating to a Toda system, the coupled ODEs can be solved exactly. In the extremal 
limit, when the asymptotic value do of a(r) is set to unity, this was done in We can 
easily redo their computations when ao is left arbitrary and express the final solutions in 
terms of two integration constants d\ and d 2 similar to the Kaluza-Klein system in the main 
body of the paper. The result can be expressed in terms of two functions 

A(r) = 12 (4 d\ - 2di d 2 ) r 4 - 12 {Ad% - d x ) r 3 + 24 d\ r 2 - 8d 2 r + 1, (B.3) 
B(r) = 6 dt r 3 - 12 d\ r 2 + 6 d 2 r - 1. (B.4) 

In terms of these functions, the solution takes the form 

exp[lo0(r)] = § ^H' a{r)= io^Qf Q r A(r)V5B (r) i/io' b ^ = ^ry ( R5 ) 

Demanding regularity in < r < 00 results in somewhat more complicated equations in 
this case, because the polynomials A(r) and B(r) are quartic and cubic and are best handled 
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Figure 4: Basin structure of this system should be contrasted to that of the KK black hole 
(figure 1) and those in [3]. 



with a computer. But since there is no rotation, this is still very much a tractable situation. 
We will not present the details, but going about the regularity conditions systematically as 
in Appendix A results in a simple result 

2d\ > c?i, dx < 0, d 2 < (B.6) 

for the attraction basin. The boundaries of attraction are given by 

2d\ = d u dx = 0, with d 2 < (B.7) 

The basin structure is presented in figure 4. 

It is possible to do a near horizon perturbation theory analogous to that done in Appendix 
D of [3] for this system as well (cf. also section 2 of the present paper), to figure out the 
trajectory of a perturbed solution along the attraction basin. The result is that the basin 
gets sliced along 

d 2 = d° 2 (B.8) 

as was also the case in the solution we discussed in section 2. Here d® is a negative constant 
that corresponds to the d 2 of the unperturbed black hole. As in the case discussed in the 
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3 L 



Figure 5: The radial profiles of the scalars in the model in Appendix 2. The strange shapes 
of the curves that start beyond the upper boundary of the basin are perfectly legitimate 
features of the plot. 



main body of the paper, the leading scalar perturbation (call it 03 ) is not at 0(r), here it is 
at 0(r 3 ). Perturbation theory tells us the relations 

dx = d° 1 + h 3 , d 2 = dl (B.9) 


upon matching the perturbed solution to the parameters of the new solution. Eliminating 
</>3 leads to (B.8). 

It is trivial to check that the warp factors defined via A = r 4 /a(r) 4 goes at the d\ = 
boundary as r 8//5 and at the 2d\ = dx as r 18 / 5 . This geometry is another version of the 
extremal subtracted (subttractor) geometry. 

For the specific case where d® = —2, we present the curves that arise for a selection of 
values of di, as we dial it from one side of the boundary 1, through boundary 1, through the 
basin, through boudnary 2, to the other side of boundary 2. The nature of the divergent 
solutions is somewhat different hera^l than in the other cases we have examined. But this 
is nothing pathological. The new feature is that the asymptotic values of the fields in this 
solution do not uniquely fix an ordered pair of real numbers (c^,^) in this case. In the 
previous cases, the only real pair (rfi,^) that gave rise to a given set of asymptotic data, 
was the non-divergent one. Here on the other hand, while one of the solutions falls inside 
the basin and gives rise to the expected curve, another one lies outside and is seen as a 

12 More precisely, the curves beyond the upper boundary have divergences, but they asymptote to constant 
scalah values at inifmity. 
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divergence before it crosses inside the basin. Of course this is not a regular solution, and as 
expected, in the near-horizon region which is responsible for the attraction mechanism, it 
lies outside the basin. 
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